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INTRODUCTION 


A  perforated  muzzle  brake  consists  simply  of  a  set  of  vents  drilled  through 
the  wall  of  a  cannon  near  the  muzzle, (see  Figure  1).  >Venting  reduces  the  axial 
thrust  produced  by  the  gas  at  the  muzzle,  thereby  effecting  a  decrease  in  weapon 
impulse.  However,  the  redirected  exhaust  increases  the  blast  levels  upstream  of 
the  muzzle.  From  a  designer's  perspective,  the  problem  is  to /choose  a  cannon- 
brake  system  which  yields  specified  values  of  muzzle  velocity  and  weapon  impulse 
but  minimizes  the  blast  increase. 
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Figure  1.  Schematic  drawing  of  a  perforated  muzzle  brake. 

The  magnitude  of  the  weapon  impulse  was  extensively  studied  by  Dillon 
(refs  1,2),  Dillon  and  Nagamatsu  (refs  3-5),  Nagamatsu  et  al .  (refs  6,7),  and 
Carofano  (ref  8).  References  1  through  5  also  contain  shadowgraph  and  free- 
field  overpressure  data  which  characterize  the  blast  field  surrounding  the  can¬ 
non. 

The  computation  of  the  blast  field  is  challenging  because  the  flow  is  tran¬ 
sient,  three-dimensional,  and  must  be  computed  to  distances  considerably  beyond 
the  breech  if  realistic  pressure  histories  near  the  breech  are  to  be  obtained 
for  reasonable  periods  of  time.  The  time  requirement  is  necessary  because  the 
peak  pressure  at  a  particular  location  is  not  always  achieved  in  the  initial 
portion  of  the  blast  wave.  To  render  the  calculation  tractable,  Carofano  (ref 
9)  exploited  a  number  of  features  of  the  flow.  A  brief  description  follows  (see 
also  References  6  and  8). 
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When  the  propellant  gas  expands  through  the  brake,  an  asymmetric  pressure 
distribution  develops  in  each  vent  with  the  highest  pressures  acting  on  the 
downstream  surface.  To  calculate  the  blast  field,  the  flow  through  each  vent  is 
required  at  each  instant  of  time  during  tube  blowdown.  Because  the  flow  is 
three-dimensional,  it  is  not  practical  to  obtain  the  complete  solution  with  a 
transient  calculation.  Fortunately,  the  flow  contains  many  features  which  per¬ 
mit  a  vigorous  simplification  of  the  problem. 

First,  because  of  the  large  volume  of  the  gun  tube,  the  blowdown  process 
takes  on  the  order  of  tens  of  milliseconds,  while  the  three-dimensional  calcula¬ 
tions  indicate  that  the  flow  in  a  vent  is  established  in  a  fraction  of  a  mi11i- 
second.  Therefore,  the  latter  can  be  treated  as  quasi-steady  and  only  the  flow 
inside  and  outside  of  the  tube  must  be  considered  as  time-dependent. 

Secondly,  in  the  applications  of  interest,  the  flow  is  supersonic  as  it 
enters  the  brake  and,  due  to  the  venting,  it  expands  to  higher  Mach  numbers  as 
it  travels  downstream.  Also,  because  of  the  high  tube  pressures,  the  gas  exits 
each  hole  at  near  sonic  or  supersonic  velocities  over  most  of  the  exit  plane 
area.  Experience  has  shown  that  the  flow  is  rather  insensitive  to  the  outflow 
boundary  condition  over  the  remaining  subsonic  portion.  Thus,  the  flow  at  a 
particular  vent  location  is  not  influenced  by  events  occurring  farther 
downstream  or  outside  of  the  tube.  It  depends  solely  on  the  conditions  in  the 
tube  upstream  of  the  vent.  In  particular,  it  was  shown  (refs  6,8,9)  that  the 
flow  is  completely  described  by  the  upstream  Mach  number,  the  specific  heat 
ratio  and  the  covolume  of  the  gas,  and  the  vent  geometry.  One  solution  with 
these  parameters  specified  is  valid  for  all  upstream  pressures  and  densities. 
Thus,  while  a  wide  range  of  physical  states  are  encountered  during  blowdown, 
only  a  few  three-dimensional  solutions  are  required  to  describe  them. 
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Data  from  these  solutions  are  used  to  obtain  average  values  of  density  and 
pressure  along  with  the  mass  and  momentum  fluxes  in  the  exit  plane  of  the  vent. 
The  averages  are  dimensionless  functions  of  the  parameters  that  appear  in  the 
three-dimensional  solutions  and  are  used  to  couple  the  interior  and  exterior 
f 1 ows . 

The  transient  flow  inside  of  the  tube  is  calculated  using  the  one¬ 
dimensional  Euler  equations  with  a  source  term  included  to  represent  the  venting 
at  the  tube  wall.  This  is  constructed  from  the  mass  flux  function  and  the  local 
conditions  prevailing  in  the  tube  at  a  given  instant. 

The  transient  flow  outside  of  the  tube  is  treated  as  axisymmetr ic .  The 
large  number  of  vents  typical  of  such  brakes  and  their  symmetrical  placement 
around  the  tube  makes  this  feasible.  Since  the  area  of  each  vent  represents 
only  a  portion  of  the  local  tube  area,  the  averaged  variables  at  the  vent  exit 
have  to  be  adjusted  to  provide  an  appropriate  boundary  condition  for  the  axisym- 
metric  equations.  A  control  volume  approach  to  achieve  this  is  described  in 
Reference  9.  The  quantities  at  the  vent  exit  are  related  to  the  interior  flow 
through  the  averaged  functions  described  above.  Because  the  vent  exit  fow  -s 
supersonic,  the  extenor  boundary  condition  is  completely  determined  by  the 
local  conditions  in  the  tube. 

The  model  (ref  9)  produced  results  which  compared  favorably  with  previously 
unpublished  shadowgraphs  of  the  near-field  obtained  by  Dillon  in  his  20-mm 
experimental  program.  A  quantitative  comparison  is  made  below  of  free-field 
overpressure  data  from  that  program  taken  at  30  calibers  from  the  muzzle. 

Some  changes  have  been  made  in  the  original  model  to  make  it  more  generally 
applicable  to  cannon  designs  of  current  interest.  For  example,  the  perfect  gas 
equation  was  replaced  by  the  Abel  equation  of  state  to  more  adequately  represent 
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the  gases  at  the  pressure  levels  which  prevail  in  large  caliber  cannon.  Also, 
the  projectile  equation  of  motion  was  added  to  more  realistically  simulate  the 
flow  discharging  to  the  environment.  Previously,  the  projectile  was  restricted 
to  move  at  a  constant  velocity. 

Finally,  a  distribution  of  vents  of  variable  diameter  and  spacing  can  now 
be  accommodated.  This  is  necessary  in  large  caliber  cannon  analysis  where  the 
diameters  of  the  vents  near  the  brake  entrance  may  have  to  be  smaller  than  those 
near  the  muzzle  to  avoid  exceeding  allowable  stress  levels.  It  was  also  reeded 
to  simulate  the  upstream  venting  scheme  discussed  below. 

INITIAL  CONDITIONS 

The  starting  configuration  is  shown  in  Figure  2.  The  initial  projectile 
Dosition  is  chosen  such  that  the  precursor  shock  is  located  just  upstream  ot  -he 
vented  region.  The  state  of  the  air  between  the  shock  and  the  projectile  s 
taken  to  be  uniform  and  is  computed  from  the  projectile  velocity  at  th-s 
instant.  The  latter  is  taken  from  the  output  of  a  standard  ballistics  solut-on. 

VENTS 


D 

hit  n  n  m 

PROPELLANT  QAS 

MOVING  AIR 

STAGNANT  AIR 

BREECH  PROJECTILE  SHOCK 


Figure  2.  Starting  configuration  showing  projert-ilo 
pushing  shock  into  the  vented  region. 

The  specification  of  the  state  of  the  propellant  gas  behind  the  projectile 
is  delayed  until  its  base  reaches  the  vented  region.  During  this  interval,  an 
analytical  representation  of  the  velocity  and  position  time  histories,  taken 
from  the  ballistics  solution,  is  used  to  advance  the  projectile  and  drive  the 
numerical  solution  downstream  of  it.  The  propellant  gas  properties  are  then 
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calculated  from  the  Pidduck-Kent  limiting  solution  for  an  Abel  gas  (ref  10) 
using  the  ilistics  data  summarized  in  Table  I  (see  References  6  and  8  for  more 
details).  In  this  manner,  essentially  all  of  the  information  generated  by  the 
ballistics  solution  relating  to  the  combustion,  friction,  and  heat  transfer 
processes  is  included  in  the  starting  data  behind  the  projectile. 

TABLE  I.  STARTING  OATA  FOR  PROPELLANT  GAS 

105-mm  __  120-mm  ! 

5.92  I  6.26  ' 

5.79  13.45  | 

1466.1  1143.3 

j 

798.0  660  1  j 

427.5  402.0  ! 

I 

475.0  485.1  ! 

i 

10.5  12.0  1 

6472.9  8749.9  1 

1.24  1.23 

I 

1050.0  1035.0  j 

4.76  3.69  ! 

i 

A  comment  is  necessary  regarding  the  starting  data  in  Table  I  for  the  20-mm 
cannon.  The  vents  were  actually  part  of  an  11  84-cm  extension  rather  than 
integral  with  the  tube.  In  the  bare  muzzle  case  then,  the  projectile  base 
travel  was  only  143.0  cm.  The  precursor  shock  was  initially  placed  at  the 
muzzle  rather  than  upstream,  as  depicted  in  Figure  2,  and  the  propellant  gas 
properties  were  specified  when  the  projectile  base  reached  the  muzzle. 


Parameter/Cannon 

20-mm 

Propellant  mass,  kg 

0.0389 

Projectile  mass,  kg 

0.0980 

Projectile  velocity,  m/sec 

1045.0 

Projectile  base  pressure,  atm 

287.0 

Projectile  base  position,  cm 

143.0 

Projectile  base  travel,  cm 

154.84 

Bore  diameter,  cm 

2.0 

Gun  chamber  volume,  cm3 

41.7 

Specific  heat  ratio 

1.25 

Covolume,  cm3 /kg 

982 . 0 

Vent  area  ratio  of  brake 

6.69 

5 


VENT  PATTERNS 


Each  brake  had  12  columns  of  vents  uniformly  spacea  around  the  tube  circum¬ 
ference  (columns  run  parallel  to  the  tube  axis).  One  column  for  each  brake  -c 
shown  schematically  in  Figure  3.  All  dimensions  are  scaled  by  the  cannon  pore 
diameter  to  facilitate  comparison.  The  vent  area  ratio,  defined  as  the  rat  o  rf 
the  total  vent  area  to  the  cannon  bora  area,  is  given  as  the  last  entry  in  'able 
I . 

0°0000000000000°  MUZZLE 
2Q-MM  CANNON  VENT  PATTERN 

OOOOOOOOOOOO  MUZZLE 

105-MM  CANNON  VENT  PQTTERN 

O  OOOOOOO  MUZZLE 

120-MM  CRNNQN  VENT  PQTTERN 

,  Figure  3.  Vent  pattern  for  each  cannon. 

In  the  20-mm  brake,  every  other  hole  was  offset  by  15  degrees  in  the  cir¬ 
cumferential  direction  to  produce  a  staggered  pattern.  The  code  considers  on’y 
the  vent  area  per  unit  length  of  tube  so  the  effect  of  staggering  cannot  be 
estimated.  Presumably,  this  arrangement  is  more  likely  to  produce  the  axisym- 
metric  flow  field  assumed  in  the  model  than  the  straight  patterns. 
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Vents  of  variable  diameter  and  spacing  were  used  in  the  120-mm  brake  to 
avoid  exceeding  allowable  stress  levels.  This  feature  is  considered  by  the 
code . 

A  more  complete  description  of  the  experimental  setup  for  the  20-mm  cannon, 
including  a  photograph  of  the  brake,  is  given  in  Reference  1.  Further  details 
of  the  105-mm  and  120-mm  tests  are  given  in  References  11  and  12,  respectively. 

CALCULATIONS 

Marten's  Total  Variation  Diminishing  scheme  (ref  13)  was  used  in  conjunc¬ 
tion  with  a  time-splitting  algorithm  (ref  14)  to  solve  the  Euler  equations.  The 
calculations  were  performed  on  a  Cray  X-MP/48  computer  using  a  single  processor. 
A  major  effort  was  made  to  exploit  the  vector  hardware  wherever  possible. 

A  uniform  grid  was  employed  over  a  rectangular  region  extending  60  calibers 
upstream  from  the  muzzle,  110  calibers  downstream,  and  70  calibers  radially  out¬ 
ward  from  the  tube  axis.  Beyond  this  region,  a  gradually  expanding  grid  was 
used  to  limit  memory  requirements  while  still  permitting  the  calculation  to  con¬ 
tinue.  Four  cells  were  used  across  the  tube  radius,  800  in  the  axial  direction 
and  350  in  the  radial  direction.  The  program  required  1.6  megawords  of  memory 
for  these  array  sizes. 

The  size  of  the  active  grid  is  determined  at  the  beginning  of  each  time 
step  to  eliminate  computation  in  the  undisturbed  environment.  Run  times  for 
each  configuration  are  given  in  Table  II.  More  time  steps  are  required  for  the 
bare  muzzle  cases  because  the  disturbance  which  propagates  upstream  takes  longer 
to  reach  the  pressure  gages  (see  next  section). 
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TABLE  II.  COMPUTATION  TIMES 


Cannon 

Time  Steps 

CPU  Minutes 

20-mm,  bare  muzzle 

1200 

38.8 

20-mm,  with  brake 

1000 

25.0 

105-mm,  bare  muzzle 

1200 

38.6 

105-mm,  with  brake 

1000 

28.5 

120-mm,  bare  muzzle 

1200 

37.8 

120-mm,  with  brake 

1000 

26.6 

PRESSURE  CONTOUR  PLOTS 

The  blast  fields  produced  by  each  cannon  after  600  time  steps  are  shown  in 
Figures  4  through  6.  Each  plot  is  scaled  by  the  respective  cannon  bore  diameter 
to  facilitate  comparison.  The  small  circles,  located  on  a  radius  30  calibers 
from  the  muzzle,  indicate  where  pressure  histories  were  stored  in  the  calcula 
tions  or  measured  in  the  experiments.  Note  in  Figure  4,  that  when  the  brake 
extension  was  added  in  the  20-mm  experiment,  the  gages  were  left  at  the  oosi- 
tions  they  occupied  in  the  bare  muzzle  case. 

The  principal  effect  of  venting  is  seen  to  be  the  generation  of  a  more  uni¬ 
form  blast  field  around  the  cannon.  The  disturbance  is  diminished  somewhat 
downstream  of  the  muzzle  and  considerably  strengthened  upstream. 

The  20-mm  cannon  had  a  significantly  higher  ratio  of  travel  length  to  bore 
diameter  than  either  large  cannon.  This  produced  a  precursor  flow  of  relatively 
long  duration  ahead  of  the  projectile.  In  Figure  4,  remnants  of  the  precursor 
shock  can  be  seen  upstream  of  the  muzzle  and  near  the  60-degree  gage  position  in 
the  brake  case.  The  precursor  shock  is  completely  overtaken  by  the  main  blast 
wave  for  both  large  cannon  because  of  their  relatively  shorter  barrels. 
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gure  6.  Pressure  contour  clots  for  the  120  mm  cannon  with 
and  without  venting. 


What  other  differences  exist  in  the  various  blast  fields  are  due  mainly  to 
the  variations  in  projectile  base  pressure  and  velocity  or  the  brake  geometry. 
Some  20-mm  experiments  are  planned  that  will  employ  brakes  which  are  geometri¬ 
cally  similar  to  those  being  used  in  ongoing  105-mm  and  120-mm  tests.  The 
question  of  scaling  the  blast  field  will  be  addressed  when  these  data  become 
available. 

OVERPRESSURE  RESULTS 

In  Figure  7,  experimental  free-field  overpressure  data  are  given  for  the 
20-mm  cannon.  Each  cluster  of  data  symbols  is  the  result  of  four  shots.  The 
zero  angle  coincides  with  the  projectile  flight  path. 

The  solid  and  dashed  lines  in  the  figure  correspond  to  calculations  made 
with  the  Abel  and  perfect  gas  equations,  respectively.  The  projectile  base 
pressure  for  this  cannon  is  less  than  300  atmospheres  so  the  covolume  term  in 
the  Abel  equation  had  only  a  modest  effect  on  the  results.  The  predictions  are 
generally  in  good  agreement  with  the  data.  The  only  exception  occurs  at  the 
10-degree  position  with  the  brake  in  place. 

The  results  in  Figure  8  are  for  the  105-mm  cannon.  In  this  case,  the  base 
pressure  is  near  800  atmospheres  and  the  value  of  using  the  Abel  equation  is 
evident.  The  predictions  lie  somewhat  above  the  data  forward  of  the  muzzle,  but 
this  is  not  a  characteristic  of  the  model.  The  comparison  in  Figure  9  for  the 
120-mm  cannon  and  a  base  pressure  of  660  atmospheres,  shows  more  satisfactory 
agreement  at  these  locations. 

Of  some  concern  is  the  tendency  of  the  model  to  underpredict  the  data  at 
the  150-degree  location  for  both  large  cannon.  The  experimental  peak  may  be  due 
to  a  wave  reflected  off  the  ground  or  the  vehicle.  The  cause  is  being  investi¬ 
gated. 
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Figure  8.  Comparison  of  model  predictions  with  experimental 

free-field  overpressure  data  for  the  105-min  cannon. 


DEGREES 


Note  that,  in  general,  the  peak  overpressure  is  somewhat  lower  near  the 
projectile  flight  path  and  rises  to  a  maximum  farthor  off  axis.  This  can  be 
explained  by  reference  to  the  contour  plots.  Near  the  axis,  the  disturbance  is 
influenced  by  the  relatively  weak  projectile  bow  wave.  Farther  off  axis,  the 
disturbance  is  due  to  the  strongest  part  of  the  main  blast  wave.  The  only 
exception  occurs  in  the  20-mm  brake  case  where  the  gage  nearest  the  axis  is 
struck  directly  by  the  main  blast  wave  because  of  its  position  relative  to  the 
exit  plane  of  the  brake  extension. 

The  overpressure  plots  for  the  large  cannon  indicate  that  venting  decreases 
the  strength  of  the  blast  wave  downstream  of  the  muzzle  and  increases  it 
upstream.  For  the  20-mm  cannon,  the  reduction  downstream  was  not  as  pronounced, 
which  is  again  associated  with  the  brake  extension.  In  a  design  situation,  the 
tube  will  have  to  be  lengthened  somewhat  to  maintain  the  desired  muzzle  veloc¬ 
ity,  but  the  addition  will  be  less  than  the  vented  length  because  the  projectile 
continues  to  accelerate  through  this  region.  The  trend  of  the  ovemressure 
measurements  will  then  lie  somewhat  intermediate  between  the  extremes  observed 
here.  In  any  event,  the  upstream  pressure  levels  will  increase.  A  method  to 
limit  the  rise  is  discussed  in  the  next  section. 

BLAST  REDUCTION  NEAR  THE  BREECH 

In  the  experiments,  all  of  the  vents  were  located  near  the  muzzle.  The 
question  arises,  "Could  another  arrangement  reduce  the  blast  levels  near  the 
breech  without  producing  significant  changes  in  weapon  impulse  or  projectile 
velocity?"  Several  patterns  have  been  considered. 

The  obvious  choice,  increasing  the  vent  spacing  to  spread  the  brake  dis¬ 
turbance  over  a  larger  area,  raised  the  pressure  levels  at  the  breech  signifi¬ 
cantly.  Displacing  so  many  of  the  vents  upstream  simply  moves  the  source  of  the 
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disturbance  closer  to  the  area  of  concern. 

The  most  successful  change  involved  moving  just  one  or  two  rows  of  vents 
upstream  while  leaving  the  remaining  vents  at  the  muzzle.  Overpressure  predic¬ 
tions  are  shown  in  Figure  10  for  the  region  aft  of  the  muzzle  of  a  105-mm  cannon 
at  a  50-caliber  radius  (a  different  round  was  used  in  these  calculations  than 
the  one  described  above).  The  bare  muzzle  data,  represented  by  the  square  sym¬ 
bols,  decrease  monotonically  from  muzzle  to  breech.  Adding  12  rows  of  vents 
near  the  muzzle  produces  the  opposite  trend,  as  indicated  by  the  circles. 
Displacing  two  of  these  rows  ten  calibers  upstream  from  the  muzzle  produces  a 
significant  reduction  in  blast  near  the  breech,  as  indicated  by  the  triangles. 

DATA  AT  50  CALIBERS  FROM  MUZZLE 

D  BARE  MUZZLE 
O  12  ROWS  AT  MUZZLE 
A  10  ROWS  AT  MUZZLE 


2  ROWS  10  CALIBERS  UPSTREAM 

9 


Figure  10.  Overpressure  predictions  aft  of  the  muzzle 
of  a  105-mm  cannon  at  a  50-caliber  radius. 
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The  scheme  appears  to  work  for  two  reasons.  First,  because  only  a  small 
number  of  vents  are  moved  upstream,  the  blast  wave  produced  by  them  is  rela¬ 
tively  weak.  Secondly,  the  flow  field  associated  with  them  interferes  with  the 
propagation  of  the  blast  wave  produced  by  the  remaining  vents.  The  result  is 
that  two  waves  arrive  at  the  breech  rather  than  one.  By  moving  the  right  number 
of  vents  the  correct  distance  from  the  muzzle,  the  strength  of  each  wave  can  be 
minimized. 

As  noted  above,  a  vented  tube  must  be  somewhat  longer  to  maintain  the 
desired  projectile  velocity.  Displacing  some  of  the  vents  upstream  adds  another 
fraction  of  a  caliber.  However,  because  the  upstream  vents  work  at  a  higher 
pressure  level,  fewer  vents  are  required  to  match  the  impulse  reduction  obtained 
with  the  unsplit  design.  In  the  example  above,  only  11  rows  are  needed  rather 
than  the  original  12. 

The  scheme  will  be  tested  in  the  laboratory  using  a  20-mm  cannon  and  in  the 
field  using  a  105-mm  cannon.  More  complete  details  of  the  calculations  will  be 
given  in  a  future  report  when  the  test  results  become  available. 

CONCLUSIONS 

The  model  predictions  are  in  satisfactory  agreement  with  available  over¬ 
pressure  data  for  small  and  large  caliber  cannon.  The  covolume  correction  in 
the  Abel  equation  significantly  improves  the  results  for  the  latter.  More  data 
will  be  available  in  the  near  future  for  further  comparison. 

Upstream  venting  shows  considerable  promise  as  a  method  of  reducing  blast 
levels  in  the  breech  area  while  maintaining  specified  values  of  projectile 
velocity  and  weapon  impulse.  Laboratory  and  field  testing  of  the  scheme  are 
planned. 
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